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ABSTRACT 

I shortly review the present status of the theoretical calculations of e'/e and the 
comparison with the present experimental results. I discuss the role of higher 
order chiral corrections and in general of non-factorizable contributions for the 
explanation of the AI — 1/2 selection rule and direct CP violation in kaon 
decays. Still lacking satisfactory lattice calculations, analytic methods and 
phenomenological approaches are helpful in understanding correlations among 
theoretical effects and experimental data. Substantial progress from lattice 
QCD is expected in the coming years. 



1 Introduction 

The results obtained in the last few years by the NA48 0) and the KTeV i) col- 
laborations have marked a great experimental achievement, establishing some 



35 years after the discovery of CP violation in the neutral kaon system U> the 
existence of a much smaller violation acting directly in the decays: 



„,„•../ , / (15.3 ± 2.6) x 10- , 
tt-eie ,e,-\ (20.7 ± 2.8) x 10~ 4 fwnar\ ! 



(NA48) 
(KTeV). 



The average of these results with the previous measurements by the NA31 
collaboration at CERN and by the E731 experiment at Fermilab gives 

Re(e'/e) = (17.2 ± 1.8) x 10~ 4 . (2) 

While the Standard Model (SM) of strong and electroweak interactions pro- 
vides an economical and elegant understanding of indirect (e) and direct (e') 
CP violation in term of a single phase, the detailed calculation of the size of 
these effects implies mastering strong interactions at a scale where perturba- 
tive methods break down. In addition, direct CP violation in K — > 7T7t decays 
arises from a detailed balance of two competing sets of contributions, which 
may hopelessly inflate the uncertainties related to the relevant hadronic matrix 
elements in the final outcome. All that makes predicting e'/e a complex and 
challenging task It . 

Just from the onset of the calculation the presence in the definition of 
e'/e , written as 

e' _ 1 [ {(™)l=2\H W \K L ) ((7T7t) I=2 \H W \K S ) \ 

of given ratios of isospin amplitudes warns us of a longstanding and still un- 
solved theoretical "problem": the explanation of the AI = 1/2 selection rule. 

The AI =1/2 selection rule in K — > tttt decays is known since 45 years H) 
and it states the experimental evidence that kaons are 400 times more likely to 
decay in the I — two-pion state than in the 1 = 2 component (w = A2/A0 ~ 
1/22). This rule is not justified by any symmetry argument and, although it 
is common understanding that its explanation must be rooted in the dynamics 
of strong interactions, there is up to date no derivation of this effect from first 
principle QCD. 

Given the possibility that common systematic uncertainties may a-priori 
affect the calculation of e'/e and the AI =1/2 rule (see for instance the present 
difficulties in calculating on the lattice the "penguin contractions" for CP vi- 
olating as well as for CP conserving amplitudes H)) a convincing calculation 



of e'/e must involve at the same time a reliable explanation of the AI =1/2 
selection rule. Both observables indicate the need of large corrections to factor- 
ization in the evaluation of the four-quark hadronic transitions. Among these 
corrections Final State Interactions (FSI) play a substantial role. However, FSI 
alone are not enough to account for the large ratio of the 1 = over 1 = 2 
amplitudes. Other sources of large non-factorizable corrections are therefore 
needed for the CP conserving amplitudes @' It, which might affect the de- 
termination of e'/e as well. As a consequence, a self-contained calculation of 
e'/e should also address the determination of the K — > hit rates. 

2 OPE: an "effective" approach 

The Operator Product Expansion (OPE) provides us with a very effective way 
to address the calculation of hadronic transitions in gauge theories. The inte- 
gration of the "heavy" gauge and matter fields allows us to write the relevant 
amplitudes in terms of the hadronic matrix elements of effective quark opera- 
tors and of the corresponding Wilson coefficients (at a scale //), which encode 
the information about those dynamical degrees of freedom which are heavier 
than the chosen renormalization scale. According to the SM flavor structure 
the AS = 1 transitions are effectively described by 

Has=i = ^V ud V; s £ [zii/i) + r W (a0] Qifr) ■ (4) 

* i 

The entries of the 3x3 Cabibbo-Kobayashi-Maskawa (CKM) matrix de- 
scribe the flavour mixing in the SM and t = —VtdY t */V u dV* s . For /x < m c 
(q = u, d, s), the relevant quark operators are: 



Qi = {s a Uf 3 ) v _ A (u l3 d a ) v _ A 
Q 2 = (sm)v-a (« d )v-A 
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Current-Current 



Gluon "penguins" (5) 



Electroweak "penguins" 



2 \°avptY-A ' 

Current-current operators are induced by tree-level W-exchange whereas the 
so-called penguin (and "box") diagrams are generated via an electroweak loop. 



Only the latter "feel" all three quark families via the virtual quark exchange 
and are therefore sensitive to the weak CP phase. Current-current operators 
control instead the CP conserving transitions. This fact suggests already that 
the connection between e'/e and the AI =1/2 rule is by no means a straight- 
forward one. 

Using the effective AS* = 1 quark Hamiltonian we can write e'/e as 



£ ' i<b G-puJ 

— = e v — — — ImAt 

e 2|e|ReA 

where 



n - - n 2 

10 



(6) 



n = ^ZiViMQih (l-fiiB) 
n 2 = ^EiViMQih > 1,1 

and (Qi) = (inr\Qi\K) . The rescattering phases 8q } 2 can be extracted from 
elastic tt-tt scattering data @) and are such that cos Jo — 0.8 and cos 82 — 1. 
Given that the phase of e (6 e ) is approximately it J A, as well as the difference 
<5o — 62, the e'/e phase <f> = ^ + 82 — Sq — 8 e turns out to be consistent with 
zero. While Gf, uj, \e\ and Rej4o are precisely determined by experimental 
data, the first source of uncertainty that we encounter in eq. ^ is the value 
of ImAt = l m (VtsVtd), the combination of CKM elements which measures CP 
violation in AS = 1 transitions. The determination of ImAt depends on B- 
physics constraints and on e & . In turn, the fit of e depends on the theoretical 
determination of Bk, the A" — K° hadronic parameter, which should be sclf- 
consistently determined within every analysis. The theoretical uncertainty on 
Bk was in the past the main component of the final uncertainty on Im A t . The 
improved determination of the unitarity triangle coming from B-factories and 
hadronic colliders has weakened and will eventually lift the dependence 
of ImAt on Bk, allowing for an experimental measurement of the latter from 
e . Within kaon physics, the decay Kl — > tt°vv gives the cleanest "theoretical" 
determination of Im A t , albeit representing a great experimental challenge. At 
present, a typical range of values for ImA t is (0.94 — 1.60) x 10~ 4 0. 

We come now to the quantities in the square brackets. While the calcu- 
lation of the Wilson coefficients is well under control, thanks primarily to the 
work done in the early nineties by the Munich HI) and Rome 111) groups, 
the evaluation of the "long-distance" factors in eq. ^ is the crucial issue for 
the ongoing calculations. The isospin breaking (IB) parameter fiiB, gives at 
the leading-order (LO) in the chiral expansion a positive correction to the A2 



amplitude (proportional to Aq via the ir° — -q mixing ) of about 0.13 0). At the 
next-to-leading order (NLOJthe full inclusion of the tt° — r\ — rf mixing lift the 
value of f^iB to 0.16 ±0.03 113). On the other hand, the complete NLO calcula- 
tion of IB effects beyond the tt° — r\ — rj mixing (of strong and electromagnetic 
origin, among which the presence of AI = 5/2 transitions) involves a number of 
unknown NLO chiral couplings and is presently quite uncertain. Dimensional 
estimates show that IB effects may be large and affect e'/e sizeably in both 
directions lil) . Although a partial cancellation of the indirect ( Alu) and direct 
Af2iB NLO isospin breaking corrections in eq. ^| may reduce their final numer- 
ical impact on e'/e , we must await for further analyses in order to confidently 
assess their relevance. At present one may use fiie = 0.10 ± 0.20 
as a conservative estimate of the IB effects. 

The final basic ingredient for the calculation of e'/e is the evaluation of 
the K — > irir hadronic matrix elements of the quark operators in eq. |[ A simple 
albeit naive approach to the problem is the Vacuum Saturation Approximation 
(VSA), which is based on two drastic assumptions: the factorization of the four 
quark operators in products of currents and densities and the saturation of the 
intermediate states by the vacuum state. As an example: 

(Tr + Tr-\Qe\K°) = 2 (it- \u l5 d\0) (tt + \su\K°) - 2 (ir + it' \dd\0) (0\s-f 5 d\K°) 
+ 2 [(0|ss|0) - (0|dd|0)] (ir+ir-\s-f 5 d\K°) (8) 

The VSA does not exhibit a consistent matching of the renormalization scale 
and scheme dependences of the Wilson coefficients and it carries potentially 
large systematic uncertainties El) . On the other hand it provides useful insights 
on the main features of the problem. A pictorial summary of the relative 
weights of the contributions of the various operators to e'/e , as obtained in 
the VSA, is shown in Fig. |]. 

As we have already mentioned, CP violation involves loop-induced op- 
erators (Q3 — Qio)- From Fig. |l| one clearly notices the potentially large 
cancellation among the strong and electroweak sectors and the leading role 
played by the gluonic penguin operator Qq and the electroweak operator Qg. 
Tipical range of values for e'/e , obtained using the VSA, are shown in Fig. || 
together with the three most updated predictions available before 1999 (when 
the first KTcV and NA48 results became known) S H> © . The fact that 
the cancellation among the strong and electroweak sectors turns out to be quite 
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Figure 1: Anatomy of e'/e in the Vacuum Saturation Approximation. In light 
(dark) gray the positive (negative) contributions of the effective four-quark 
operators are shown with proportional weight. 

effective (in the VSA) warns us about the possibility that the uncertainties in 
the determination of the relevant hadronic matrix elements may be largely am- 
plified in the calculation of e'/e . It is therefore important to asses carefully 
the approximations related to the various parts of the calculations. In par- 
ticular, the analysis of the problem suggests that factorization may be highly 
unreliable. 

3 Beyond Factorization 

The dark gray bars in Fig. || depict the results of three calculations of e'/e which 
are representative of approaches that (in principle) allow us to go beyond 
naive factorization. They are based from left to right on the large N c ex- 
pansion EH EH), on lattice regularization EH' EH\ and on phenomenological 
modelling of low-energy QCD (the chiral quark model) c3> c3 EH) . 

The experimental and theoretical scenarios have changed substantially 
after the first KTeV data and the subsequent NA48 results. Fig. || shows 
the present experimental world average for e'/e compared with the revised 
or new theoretical calculations that appeared during the last year. Without 



all represent attempts to incorporate non-perturbative information into the 



entering into the details of the results (for a short summary see 
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PRE-DICTIONS 

Figure 2: The 1-a results of the NA31 and E731 Collaborations (early 90's) are 
shown by the gray horizontal bands. The old Miinchen, Roma and Trieste the- 
oretical predictions for e'/e are depicted by the vertical bars with their central 
values. For comparison, the VSA estimate is shown using two renormalization 
schemes. 



calculation of the hadronic matrix elements, whether their are based on the 
large N c expansion (Miinchen HI), Dortmund H^), Beijing HI), Taipei H), 



Valencia 
Trieste [ 



phenomenological modelling of low-energy QCD (Dubna 



Lund 



CD Sum Rules (Montpcllicr Eli)) or, finally, on 
lattice regularization (Roma d), CP-PACS H), RBC EH). 

Overall most of the theoretical calculations are consistent with a non- 
vanishing positive effect in the SM (with the exception of the recent lattice 
results on which I will comment shortly). 

At a closer look however, if we focus our attention on the central values, 
many of the predictions prefer the 10 -4 regime, whereas only a few of them 
stand above 10~ 3 . Is this just "noise" in the theoretical calculations? Without 
entering the many details on which the estimates are based, most of the afore- 
mentioned difference can be explained in terms of a single effect: the different 
size of the hadronic matrix element of the gluonic penguin Qq as obtained in 
the various approaches. In turn, this can be understood in terms of sizeable 
higher order chiral contributions (NLO in the 1/N C expansion) to the 1 = 
amplitudes. 
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Figure 3: Recent theoretical calculations of e'/e are compared with the 
combined \-a average of the NA31, E731, KTeV and NA48 results (e'/e = 
17.2 ± 1.8 x 10~ 4 ), depicted by the gray horizontal band. 



This effect was stigmatized well before the latest experimental round by 
the work of the Trieste group If! 



ID , and appears clearly in the comparison 
of the leading l/N c and lattice results with the chiral quark model analysis 
in Fig. ^. The chiral quark model approach, together with the fit of the CP 
conserving amplitudes which normalizes phenomenologically the matching and 
the model parameters, allows us to carry the calculation of the hadronic ma- 
trix elements beyond the leading order in the chiral expansion (including the 
needed local counterterms). Non-factorizable chiral contributions (missing in 
the leading l/N c or lattice calculations) were shown to produce a substantial 
enhancement of the I — transitions thus lifting the expectation of e'/e at the 
1(T 3 level. 

Since then a number of groups have attempted to improve the calculation 
of K — > hit matrix elements in a model independent way. Table [3] presents a 
comparison of different calculations of the relevant matrix elements. Due to 
the leading role played by Q§ and Q% we may write a simplified version of eq. |6[ 
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which although "not be used for any serious analysis" gives an effective 
and practical way to test and compare different calculations. 



The B-factors Bi = (Qi) / (Qj) V sA represent a convenient parametriza- 
tion of the hadronic matrix elements, albeit tricky, in that their values are in 
general scale and renormalization-scheme dependent, and a spurious depen- 
dence on the quark masses is introduced in the result whenever quark densities 
are involved. The latter is the case for the Qq and Qs penguins. As a conse- 
quence the VSA normalization may vary from author to author thus introduc- 
ing systematic ambiguities. By taking the VSA matrix elements at the scale 
(j, = 2 GeV we obtain @) 

<(^) 2 |Q8|^ ) VSA = V6 f m* K (m s + m d )- 2 ~ 1.1 GeV 3 , 

(HQ 6 |K°) VSA /((™) 2 |Q 8 |!^) VSA = -2\/2 {fx - U)/U =s -0.63 , 10 

where I have used (m s + md)(2 GeV) = 110 MeV and the chiral value / = 86 

MeV for the octet decay constant. 

(2) 

It is known that Bq and B 8 are perturbatively very weakly dependent 
on the renormalization scale 0). Therefore it makes sense to compare the B's 
obtained in different approaches, where the matrix elements (Qi) are computed 
at different scales. The results for the relevant penguin matrix elements coming 
from various approaches are collected in Table |3|, paying care to normalizing 
the data in a homogeneous way (as far as detailed information on definitions 
and renormalization schemes was available) . 

As a guiding information, taking ImAt = 1-3 x 10~ 4 , the present experi- 
mental central value of e'/e is reproduced by Bg (1 — flm) ~ 0.4 Bg w 1. 

The most important fact is the first evidence of a signal in lattice calcula- 
tions of (tt\Q 6 \ K), obtained by the CP-PACS EH and RBC EH collaborations. 
Both groups use the Domain Wall Fermion approach which allows to control 
the chiral symmetry on the lattice as a volume effect in a fifth dimension. This 
approach softens in principle the problem of large power subtractions which 
affects the lattice extraction of I = amplitudes (penguin contractions). Still 
only the (7r|Qj|iif) transition is computed on the lattice and LO chiral perturba- 
tion theory is used to extrapolate it to the physical amplitude. The two groups 
obtain comparable values of the Qe and Qs matrix elements leading both to 
a negative e'/e (and do not agree on the CP conserving 1 = amplitude). 
On the other hand the calculations are at an early stage and do not include 
higher order chiral dynamics which may be responsible for the enhancement of 
1 = amplitudes (as large N c approaches beyond LO and the Chiral Quark 
Model strongly suggest). The SPQcdR collaboration has reported a result for 



Table 1: Comparison of various calculations of penguin matrix elements. The 
data marked by the star are rescaled by a factor y/3/2, to account for a different 
definition of the isospin matrix elements. 



Method 


Br fNDRl 


R (2) (NDRl 




ijtXh Lice yu \\ r 7 / \ * 7T j 


\ u.o 




pp pApc |38|) 


Lattice [UWt, K — > it) 


~ U.4 


~ 1 


nnp 59b 


Lattice (if ^ tt + X PT) 




0.58 ±0.06 * 


APE© 


Lattice (if -> tt + X PT) 




0.56 ±0.07 * 


SPQcdR© 
SPQcdR© 


Lattice (K — > 7T7r) 




0.64 ±0.07 * 


Large iV c +LMD (%-limit) 




2.6 ±0.8 * 


Marseille H) 


Dispersive+data (x-limit) 




2.5 ±0.8 * 


Amherst El) 


Dispersive+data (x-limit) 




1.4 ±0.6 


Lund© 


Large N c + data 


1.0 ±0.3 


0.8 ±0.2 


Munich H) 


NLO 1/N e CHPT 


1.5 - 1.7 


0.4-0.7 


Dortmund © 


NLO 1/N C ENJL (x-limit) 


2.9 ±0.5 


1.5 ±0.2 


Lund © 


NLO xQM + xPT 


1.5 ±0.4 


0.84 ±0.04 


Trieste 13) 


Large 7V C + FSI 


1.55 ±0.10 


0.92 ±0.03 


Valencia © 



the Qs matrix element from direct calculation of the K — > 7T7t amplitude on 
the lattice. This result agrees with previous lattice data, albeit it does not yet 
include the chiral corrections relevant to quenching and to the extrapolation to 
the physical pion mass El) . 

Among the analytic approaches important results have been obtained 
using data on spectral functions in connection with QCD sum rules and dis- 
persive relations in the attempt to obtain model independent information on 
the relevant matrix elements. These approaches have produced as of today 
calculations of Q% (in the chiral limit) which are subtantially larger than the 
factorization (and lattice) results. While there is still disagreement among the 
different analysis, we must await the calculation of the Qq matrix element and 
a quantitative assessment of chiral breaking effects before drawing conclusions 
on these as well as lattice results. 

Calculations which sofar have allowed for the determination of all relevant 
parameters, based on chiral perturbation theory and/or models of low-energy 
QCD, have shown the crucial role of higher order non-factorizable corrections 



in the enhancement of the 1 = matrix elements H' —< HI). Chiral loop 
corrections drive the final value of e'/e in the ballpark of the present data. How- 
ever the calculation of higher order chiral effects cannot be fully accomplished 
in a model independent way due to the many unknown NLO local couplings. 
In the chiral quark model approach all needed local interactions are computed 
in terms of quark masses, meson decay constants and a few non-perturbative 
parameters as quark and gluon condensates. The latter are determined self- 
consistently in a phenomenological way via the fit of the CP conservin g K — > irn 
amplitudes, thus encoding the AI =1/2 rule in the calculation Sll^). The 
analysis shows that the role of local counterterms is subleading to the chiral logs 
when using the Modified Minimal Subtraction (as opposed to the commonly 
used Gasser-Leutwyler prescription). The phenomenological fit is crucial in 
stabilizing the numerical prediction IH) . The fact that the model parame- 
ters (quark and gluon condesates, constituent quark mass) turn out to be in 
the expected range, shows that the explicitly included chiral (and 1/N gluon 
condensate) corrections represent the largest non-factorizable effect. 

Among higher order corrections FSI play a leading role. As a matter of 
fact, one should in general expect an enhancement of e'/e with respect to the 
naive VSA due to FSI. As Fermi first argued El), in potential scattering the 
isospin 1 = two-body states feel an attractive interaction, of a sign opposite 
to that of the 1 = 2 components thus affecting the size of the corresponding 
amplitudes. This feature is at the root of the enhancement of the 1 = ampli- 
tude over the 1 = 2 one and of the corresponding enhancement of e'/e beyond 
factorization. An attempt to resum these effects in a model independent way 
has been worked out by the authors of ref. E^, using a dispersive approach 
a la Omnes-Mushkelishvili El' E3). Their analysis shows that resummation 
does not substantially modify the one-loop perturbative result and, as it ap- 
pears from Table |3|, a 50% enhancement of the gluonic penguin matrix element 
is found over the factorized result. However, the calculation suffers from a 
sistematic uncertainty due to the indetermination of the off-shell amplitude 
which is identified with the large iV c result El). Even when the authors in 
the most recent work match the dispersive resummation with the on-shell per- 
turbative one-loop calculation, thus including 1/N C effects, again a systematic 
uncertainty remains in the unknown polinomial parts of the local chiral coun- 
terterms. Therefore, a model-independent complete calculation of chiral loops 



for K — ► 7T7T is still missing. 

Finally, it has been recently emphasized El) that cut-off based approaches 
should pay attention to higher-dimension operators which become relevant for 
matching scales below 2 GeV and may represent one of the largest sources of 
uncertainty in present calculations. The results of refs. El' El) include these 
effects. The calculations based on dimensional regularization may be safe if 
phcnomenological input is used in order to encode in the relevant hadronic 
matrix elements the physics at all scales (this is done in the Trieste approach) . 

In summary, while model dependent calculations suggest no conflict be- 
tween theory and experiment for e'/e , a precise and "pristine" prediction of 
the observable is still quite ahead of us. 



4 Outlook and Conclusions 



Higher-order chiral corrections are taking the stage of K — > ttt: physics. They 
are needed in order to asses the size of crucial parameters (as fiis) and the 
effect of non-factorizable contributions in the penguin matrix elements. 

Lattice, as a regularization of QCD, is the first-principle approach to the 
problem. However, lattice calculations still heavily depend on chiral perturba- 
tion theory HI) . Presently, very promising developments are being undertaken 
to circumvemt the technical and conceptual shortcomings related to the calcu- 



lation of weak matrix elements llif Ell). Among those are the Domain Wall 
Fermion approach EH) which allows us to decouple the chiral symmetry from 
the continuum limit, and the very interesting observation that the Maiani- 
Testa theorem HI) can be overcomed using the fact that lattice calculations 
are performed in finite volume @) , thus allowing for the direct calculation of 
the physical K — ► tttt amplitude on the lattice. All these developments need a 
tremendous effort in machine power and in devising faster algorithms. Prelim- 
inary results for lattice calculations of both e'/e and the AI =1/2 selection 
rule are already available and others are currently under way 

In the meantime analytical and semi-phenomenological approaches have 
been crucially helpful in driving the attention of the community on some sys- 
tematic short-comings of "first-principle" calculations. The amount of theo- 
retical work triggered by the NA48 and KTeV data promises rewarding and 
perhaps exciting results in the forthcoming years. 
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